Exact solution to the averaging problem in cosmology 
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The exact solution of a two-scale Buchert average of the Einstein equations is derived for an 
inhomogeneous universe which represents a close approximation to the observed universe. The two 
scales represent voids, and the bubble walls surrounding them within which clusters of galaxies are 
located. As described elsewhere [New J. Phys. 9 (2007) 377], apparent cosmic acceleration can be 
recognised as a consequence of quasilocal gravitational energy gradients between observers in bound 
systems and the volume average position in freely expanding space. With this interpretation, the 
new solution presented here replaces the Friedmann solutions, in representing the average evolution 
of a matter-dominated universe without exotic dark energy, while being observationally viable. 

PACS numbers: 98.80.-k, 98.80.Es, 98.80.Jk, 95.36.+X 

At the time of last-scattering the distribution of mat- 
ter in the universe was very smooth, given the evidence of 
the cosmic microwave background (CMB). At the present 
epoch, by contrast, the universe is very lumpy on scales 
less than 100-300Mpc, with clusters of galaxies strung in 
filaments and bubbles surrounding huge voids. Some 40- 
50% of the present volume Jj] of the universe is in voids 
of order 30/i~^Mpc in diameter, h being the dimension- 
less Hubble parameter, = 100/ikm sec~^ Mpc^^, and 
when smaller and larger voids (2] are taken into account, 
then the observable universe is "void-dominated" . 

In spite of present-day inhomogeneity, a broadly 
isotropic Hubble flow is observed when one averages on 
sufficiently large scales. This is taken as justification for 
assuming that cosmic evolution can be modelled by the 
Friedmann equation for a smooth fluid, despite the obvi- 
ous observational evidence that galaxies are not smoothly 
distributed. To achieve agreement with a number of inde- 
pendent observations within the smooth fluid paradigm, 
dark energy has been included in the standard cosmolog- 
ical model, posing a foundational mystery for physics. 

In recent years, a number of cosmologists have ques- 
tioned whether the observations, which have been inter- 
preted as cosmic acceleration, might in fact be accounted 
for by taking more care in deriving the geometry that 
comes from averaging the actual inhomogeneous matter 
distribution. In particular, the geometry which arises 
from the time evolution of an initial average of the mat- 
ter distribution does not generally coincide, at a later 
time, with the average geometry of the full inhomoge- 
neous matter distribution evolved via Einstein's equa- 
tions Whether or not the resulting "back-reaction" 
of inhomogeneities on the average geometry can be large 
enough to explain effects usually attributed to cosmic 
acceleration from dark energy in Friedmann-Lemaitre- 
Robertson- Walker (FLRW) models, has been the subject 
of intense debate. (See 01 and references therein.) 

In this Letter I will derive the general exact solution 
to the Buchert equations Q for the two-scale model in- 
troduced in ref . [5| , yielding a simple new observationally 



viable model of the universe. The observational claim is 
based on my proposal that the debate about dark energy 
from structure formation ^ can be resolved by careful 
consideration of the operational interpretation of mea- 
surements in cosmology from first principles [Hj- This is 
necessary when averaging an inhomogeneous cosmology, 
since in general the rods and clocks of observers will be 
calibrated differently from those at an average location. 
In writing down average parameters one must define how 
they are related to our measurements operationally. As- 
suming the Copernican principle, the fact that we observe 
an almost isotropic CMB means that other observers 
should also measure an almost isotropic CMB. However, 
it does not demand that such observers measure the same 
mean CMB temperature as we, nor the same angular 
scale for the Doppler peaks in the anisotropy spectrum. 
Significant differences can arise due to gradients in spa- 
tial curvature and associated gravitational energy. 

In general relativity space is dynamical and can carry 
energy and momentum. By the strong equivalence prin- 
ciple, since the laws of physics must coincide with those of 
special relativity at a point, it is only internal energy that 
can be localised in an energy-momentum tensor on the 
r.h.s. of the Einstein equations. Thus the uniquely rel- 
ativistic aspects of gravitational energy associated with 
gradients in spatial curvature, and gradients in the ki- 
netic energy of spatial expansion, cannot be included in 
the energy momentum tensor, but are at best described 
by a quasilocal formulation. (For a review, see Q.) 

In ref. Q I propose a quantitative solution to the prob- 
lem of apparent cosmic acceleration through the techni- 
cal definition of a finite infinity scale, realising a quali- 
tative suggestion of Ellis [7| . Finite infinity replaces the 
usual notion of spatial infinity in exact asymptotically 
flat spacetimes, as the fiducial reference point for quasi- 
local gravitational energy with respect to observers in 
virialised bound systems. A universal definition of this 
scale is possible, since the initial expansion rate of the 
universe was extremely smooth at last scattering, leading 
to a true critical density, pcr, as a demarcation between 
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potentially bound and unbound systems. Due to back- 
reaction, Per, does not evolve by the Friedmann equation. 

While it has long been understood that averaging 
an inhomogeneous universe entails the dressing of av- 
erage cosmological parameters [§] through volume fac- 
tors which relate to differences in spatial curvature, the 
proposal of ref. 0] recognises that clock rates can also 
vary systematically between observers in bound systems 
within finite infinity, and a volume-average position in 
freely expanding space, due to differences in gravita- 
tional energy. By this means an iinplicit solution of the 
Sandage-de Vaucouleurs paradox [5| is possible: the lo- 
cally defined or bare Hubble parameter, H , can be uni- 
form even though voids appear to expand faster than the 
bubble walls which surround them, since cosmic clocks 
within voids tick faster on account of gravitational energy 
differences. Since our cosmological observations involve 
photons exchanged with objects in bound systems, we 
do not observe clocks in freely expanding space directly. 
Nonetheless, an ideal comoving observer within a void 
would measure a somewhat older age of the universe, and 
an isotropic CMB with a lower mean temperature and an 
angular anisotropy scale shifted to smaller angles. 

Buchert's scheme is somewhat heuristic, since it does 
not average all of the Einstein equations and requires 
an extra integrability condition to ensure closure. How- 
ever, starting from a fully covariant averaging scheme 
0, with reasonable cosmological assumptions, the corre- 
lation tensor takes the form of a spatial curvature [lo| . 
and Buchert's scheme can be realised as a consistent limit 
[llj. Following ref. the Buchert average constructed 
here is based on the two scales most relevant to the ob- 
served universe: (i) the voids which dominate the uni- 
verse at the present epoch; and (ii) finite infinity regions 
containing galaxy clusters within the filaments and bub- 
ble walls that surround voids. The local average geome- 
try at the boundary of a finite infinity region is assumed 
to be spatially flat, with the metric 



may be rewritten 



^4. 



(1) 



Within voids the metric is not given by ([T]) but is neg- 
atively curved, with local scale factor Ov • We average 
over the entire present epoch horizon volume, V = Via'^, 
where d^ = fvia^ + fwia^] fvi and /wi = l-/vi being the 
respective initial void and wall volume fractions at last 
scattering, to construct the Buchert average geometry 



ds^ = -dt^ + d^{t) df + A{f], t) dn"^ 



(2) 



Here the area function A is defined by a horizon-volume 
average 0] . The time-parameter t differs from the wall- 
time T of ll]) by the mean lapse function dt = j(t) dr. 
The geometry ([2]) is not locally isometric to the local 
geometry in either the walls or void centres. 

When the geometry ([1]) is related to the average geom- 
etry ^ by conformal matching of radial null geodesies it 
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dr' + ^[df+rl{f],T)dn' 



r 



(3) 



where = ^ — fv^^^ fwi~^^^r]w{f},T). Two sets of 
cosmological parameters are relevant: those relative to 
an ideal observer at the volume-average position in freely 
expanding space using the metric ([2]), and conventional 
dressed parameters using the metric ([3]). The conven- 
tional metric ([3]) arises in our attempt to fit a single global 
metric ([T|) to the universe with the assumption that av- 
erage spatial curvature and local clock rates everywhere 
are identical to our own, which is no longer true. 

The volume-average matter, curvature and kine- 
matic back-reaction parameters are given by fi^^ = 

87TGp^^j^dl/[3H\% = -fcv/vi'/Vv'/V[a'^'], and 

- -2 - 2 

^Iq — — /v /[9/v(l — fv)H ] respectively, where the aver- 
age curvature is due to the voids only, which are assumed 
to have fcv < 0, an overdot denotes a derivative w.r.t. 
volume-average time, t, and H = a/a is the volume- 
average or bare Hubble parameter. It satisfies 

H = fvHv + fwHvj, (4) 

where = d^/ and = a'w/ aw are the regional av- 
erage expansion rates of voids and walls, as measured by 
volume-average clocks. We define hrU) = H^/H^ < 1. 
The independent Buchert equations [3| , including the in- 
tegrability condition that ensures their closure, are 

n^.^ + + - 1, (5) 

a ^dt (f^giJ^a^) + d-^dt (n^H^d^^ = . (6) 

Conventional dressed parameters defined with respect 
to the geometry ([3]) , relevant to "wall observers" such as 
ourselves, do not satisfy a simple relation analogous to 
([5|). The conventional matter density parameter, f^^,^, 
is expected to take numerical values similar to those 
we infer in FLRW models. It differs from the bare 
volume-average density parameter, f^^j, according to 
flj^j = ^^rij^j. The mean lapse function is given by 

J=l + K'il^hr).U, (7) 

as a result of the requirement that the bare, or "locally" 
measured, expansion rate is uniform. The dressed Hub- 
ble parameter that we measure as wall observers, the 
global average over both walls and voids, is not H, but 
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The Buchert equations ([51 
pair of first order equations 



^H-T'hl- (8) 
, ^ may be reduced to the 



(1 - /v) . - i/v = V^MO^O (1 - (1 - /v) i , (9) 
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where <g; 1 is an integration constant, obtained from a 
first integral Q of eqs. ^ and (O, namely 

(l-eOf^^M(l-/v)~' = l- 
Eqs. (O and (jTU]) are readily integrated. Firstly we mul- 
tiply dSl) by H^^ (1 - /v)~^^^ a to obtain 

"(1 - f.f' 



1 d 
^dt 



^ (l-/v)' 



/3; 



with the solution 
(l-/v)^/'a-c 



MO 
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2/3 



(11) 



where a constant of integration corresponding to the ori- 
gin of time has been set to zero without loss of generality. 



Since fy 
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(1 — fv)^^^ a, we see that 
9 f .-In ■ = 

4 J Wl V 



- n f2/3 

Thus the 



where Owo = o 
local expansion rate within the wall regions is exactly 
that of an Einstein-de Sitter universe, but in volume av- 
erage time not locally measured wall time. 
Finally, we multiply pn]) by a^^H^ fv^^^a to obtain 
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where u = fv^^^a/a^ — fvi'^^'^av/a^, is proportional to 

and Ce = Ci^j^^^f^Q /flto is a constant which can ei- 
ther be positive, zero or negative depending on the initial 
value Ci. Integrating (|12p we find 
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where a 



fie) 

(13) 

and te is a constant which 
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cannot be chosen to be zero without loss of generality, as 
the time origin was already fixed in determining (jlip . 

Eqs. (|lip and (|13p constitute the general exact solution 
to the two-scale Buchert equations ([5]), ((61), regardless of 
the physical interpretation of observable quantities. In 
particular, the Buchert equations have been studied by 
a number of authors without mention as to what the 
physical interpretation of the time parameter, t, is. Here 
we will pursue the interpretational framework of ref. Q . 

A number of quantities of interest may be found di- 
rectly. Since = d^a^= 2/(3i), and 



it follows that 
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flkOU + ilj^^fl^^Ct ) (1 - /v) 



(14) 



Evaluating [([9]) -t- f„ (|TO)) ] at the present epoch, we obtain 
the following constraint on parameters 



e.)f^A/o(l - Ao) + V i^ko + f^Mo^O/vO = 1 . 

(15) 

This reduces the number of free parameters by one; we 
may take Clko as dependent, for example. The value of 
ie is also determined in terms of the other parameters by 
evaluating (|13p at the present epoch, to give 
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,(16) 



where the age of the universe in volume-average time is 
t^ = ^J^-I^, (17) 



3H. V (1 - eO^! 



A/0 

on account of pT|) . We observe that the physical interpre- 
tation of the void volume-fraction ceases to be physical 
meaningful in the limit t ^ if ie > 0, as is generally 
the case. Radiation must be included to describe the uni- 
verse at early times, and has been omitted here. Here we 
restrict attention to the matter-dominated epoch. 

The general solution is specified by four independent 
parameters, H^, ^f^^ and /vq. However, two of these 
may be further eliminated by taking priors [l^ at the 
surface of last scattering consistent with the CMB. 

Since eq. (|13p is a transcedental equation, the combi- 
nation of eqs. pT|) and ([TH]) only define a{t) and /v(i) 
implicitly in the general case. Nonetheless, at late times 
when both t and u are large, all general solutions to (|13p 
tend to the particular solution with = and te = 0. 
This particular solution is in fact a late-time tracker so- 
lution, an attractor which is insensitive to h^-i and /vi, as 
long as the observable universe is void-dominated. 

As Ce = when = 0, eq. yields a partic- 

ularly simple form for the late-time tracker solution. 



Since fvi'^^a^ 



/v 



1/3, 



OpU, we see that it corre- 



sponds to the case in which the void regions expand 
exactly as a Milne universe in volume-average time, 

flv = Ovoi, where Cvo = OL]!Qa^HJ~^'^f^i~^/'^. It fol- 
lows that hr = H^/Hy = 2/3 is a constant for this spe- 
cial solution. If we combine ([T5|) with (fTl]) evaluated 
at the present epoch, we see that only one of the pa- 
rameters ^j^jQ, i^fco and /vo is independent, and in fact, 

^Mo = 4(1 - /vo)/(2 + /vo)', ^ko = 9/vo/(2 + Uf. 
Consequently, the volume-average scale factor is 

,2/3 



ag(3i7pt) 



2 + /v. 



3/voi?o*+(l-/vo)(2 + /vo) 



1/3 



(18) 
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for the late-time tracker, while its void fraction is 



3/voi?o^+(l-/vo)(2 + /vo) 



(19) 



and the wall fraction is easily deduced from /„ = 1 — /v 
In terms of fv(t)i the mean lapse function, bare Hubble 
parameter and bare matter density have the simple forms 
7(t) = ltH{t) = ^ndn^jit) = 4(l-/v)/(2+/v)2. 

The tracker-solution dressed Hubble parameter ([5]) is 



-I 



/v(4/v + 1) 

6t 



(20) 



Other dressed parameters [5j using the metric ^ rel- 
evant to galactic wall obsevers include the redshift, z, 
and luminosity distance, — l^^a^^l + z)rv,, where 

= 7 (1 - Uf' dt'/m')il - Mt')r/^-a{t')]. For 
the tracker solution these respectively satisfy 



z + 1 
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(21) 



(22) 



This last integral can be given in a simple closed analytic 
form. These expressions are given in terms of volume- 
average time. To convert to wall time relevant to galactic 
observers, we perform the integral, t = Jq di7~^, giving 



AO. 



MO 



In 1 



(23) 



where fl 



MO 



^(1 — /vo)(2 + /vo) is the present epoch 
dressed matter density in the case of the tracker solution. 
In (f20 |) - ((22|) the parameter t should be considered to be 
implicitly defined by in terms of wall time, r. At 
late times as f — > oo, r ~ so that H ~ ~ t~^, as 
for an empty Milne universe. 

The age of the universe in volume-average time is = 
(2 + /vo)/(3ff„), and 



2(2 + /v 



(1 - /vo) 

3/vo 



In 



2(1 - /vo 



(24) 

in wall time. The dressed Hubble constant we measure, 
is related to the bare Hubble constant, H^^, by 



(4/v'o + /vo 



2(2 + /vo) 



(25) 



Finally, the volume-average tracker-solution decelera- 
tion parameter is q = ^rtj^j+2i}Q = 2 (1 — /v)^ /(2-t-/v)^, 
which begins close to the Einstein-de Sitter value, q ~ 5, 
when /v is small, and approaches g — > 0^ at late times, 
but remains positive at all times. Thus a volume-average 



observer in freely expanding space detects no cosmic ac- 
celeration. Nonetheless, a bound system observer mea- 
sures an effective dressed deceleration parameter 



(1 - /v) (8/v' + 39/v' - 12/v 
(4 + /v + 4/v^)' 



(26) 



which also begins at the Einstein-de Sitter value, 9^5, 
for small /v but then changes sign at epoch when /v ~ 
0.5867, at a zero of the cubic in ([26|) . Apparent ac- 
celeration reaches a maximum when /v ~ 0.7736 when 
q ~ —0.043, and then g ^ 0^ at late times. Thus a 
wall observer registers a late-time evolution which again 
is close to that of a Milne universe, but this time with 
apparent acceleration. 

While the initial conditions of the CMB require that 
/ir 1 at last scattering, for observationally relevant ini- 
tial conditions, the general solution pT|) . ((T3|) approaches 
the tracker solution to within 1% by a redshift of z ~ 37. 
Thus the tracker solution can be used as a very reliable 
approximation all the way back to the epoch of reioniza- 
tion. It is effectively the simplest viable generalisation of 
the Einstein-de Sitter and Friedmann models, which in- 
corporates back-reaction, and is potentially of great util- 
ity. Detailed cosmological parameter fits are presented 
elsewhere [12 L and demonstrate that the present fractal 
bubble model [5] is at the very least a serious contender for 
quantitatively solving the problem of dark energy purely 
within general relativity. The analytic solution presented 
here, and its simple tracker limit, should therefore pro- 
vide the basis for many future cosmological tests. 
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